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Two-dimensional spinless Falicov-Kimball model (FKM) with correlated hopping is studied per- 
turbatively in the limit of large on-site Coulomb interaction U. In the neutral case the effective 
Ifamiltonian in spin variables is derived up to terms proportional to U~^. Unlike the simplest FKM 
case, it contains odd parity terms (resulting from the correlated hopping) in addition to even parity 
ones. The ground-state phase diagram of the effective Hamiltonian is examined in the {a/t, h) plane, 
where a/t is a parameter characterizing strength of the correlated hopping and /i is a difference of 
chemical potentials of two sorts of particles present in the system. It appears to be asymmetric with 
respect to the change h — > —h and a new ordered phase is found for a certain interval of a/t. 



Introduction. One of the most fascinating, but still 
mysterious phenomenon observed in some materials is a 
charge ordering ||l| . Despite of its oddity, originated from 
the quantum nature of interacting electrons, it seems to 
be very likely that an inhomogeneous charge distribu- 
tion is very common effect present in strongly correlated 
electron systems Q. Then, presumably, it must be re- 
lated to a number of various phenomena found in the 
systems (such as metal-insulator phase transition, high 
Tc supecoductivity, giant magnetoresistivity, just to men- 
tion a few). The above arguments point out how much 
important is understanding of a nature of the effect and 
justify growing interest in its theoretical description 
In particular it is important to determine factors deciding 
about what sorts of charge superstructures are formed. 

One of the simplest models suitable to decribe charge 
ordered phases on a microscopic level is the Falicov- 
Kimball model, previously applied to account for the 
metal- insulator transition [||, mixed valence phenomena 
crystallization and alloy formation ||^ etc. Indeed, 
it was shown that ground state phase diagrams of the 
simplest version of the FKM have extremely rich struc- 
ture [^-|8j . The great advantage of the model is that it is 
amenable to rigorous analysis |p|,|fO|. 

However, the simplest version of the FKM, although 
non-trivial, is not able to account for all aspects of real 
experiments. This is way an objective of our studies 
is gradual inclusion of those terms that were ignored in 
the simplest version of the FKM, yet keeping the model 
tractable. 

In this contribution we investigate the FKM with so- 
called correlated hopping term added (chFKM). This 
term was already mentioned by Hubbard in ]rf[ |. More 
than a decade ago Hirsch pointed out that the term may 
be relevant in explanation of superconducting properties 
of strongly correlated electron systems (he named it the 
bond-charge interaction) [ p^ . For a last few years some 
other authors examined a role of the correlated hopping 



in the FKM ||13|, and the Hubbard model |fj,|15[l, mainly 
in a context of the metal-insulator phase transition. In- 
clusion of this term makes electron hopping rate depen- 
dent on occupation numbers of those sites between which 
an electron hops. 

One of the most difficult problems that one encounters, 
when trying to describe correlated electron systems, is a 
choice of reliable method, that enables to treat the model 
under consideration in a controllable way. Here we use 
a perturbative method valid in the large U limit, that 
permits to transform an initial Hamiltonian, having a 
small quantum part, into an effective classical one. The 
method has been reported in a series of papers by Datta, 
Fernandez, Frohlich and Rey-Bellet [l^-[l9|]. One can use 
this method to generate a perturbative series up to an ar- 
bitrarily high order in 1/U, to establish the convergence 
of the whole procedure, and - in some cases - to obtain 
phase diagrams in low (but nonzero) temperatures. This 
can be done by extending the techniques of Pirogov-Sinai 
theory pO[ | to quantum models. 

The aim of our paper is to examine properties of the 
chFKM in the perturbative regime, i.e. in the range of 
parameters where all kinds of hopping terms are small in 
comparison with the on-site Coulomb interaction term 
U . We are particularly interested in examination of how 
correlated hopping term influences charge ordering. In 
the first part of our work we derived an effective Hamil- 
tonian, that is legitimated in the large U limit. Then 
we found its ground state properties by means of the 
method of restricted phase diagrams, used previously to 
the simplest version of the FKM |@,§ . 

The model. We are dealing with two types of parti- 
cles defined on a d-dimensional simple cubic lattice Z'^: 
immobile ions and itinerant, spinless electrons. (Other 
interpretations of the model have also been considered 

1,01). 

The Hamiltonian defined on a finite subset A of Z'^ has 
the form 



f 



Ha = Ho,a + Va, 



(1) 



where 



Hq^a WxTLx - fii ^ - A^e ^ rix, (2) 



Va = - 



<xy> 



[t + a{wx + Wy)]iclcy + c)yCx) 



(3) 



Here cj, and Cx are creation and annihilation operators of 
an electron at lattice site a; € A, satisfying ordinary an- 
ticommutation relations and the corresponding number 
particle operator is Hx = c\.Cx- Wx is a classical variable 
taking values or 1. It measures the number of ions at 
lattice site x. The chemical potentials of the ions and 
electrons are and /ie, respectively, t is the electron 
hopping amplitude between empty sites and a is the cor- 
related hopping constant. The symbol < xy > denotes 
an orderless pair of nearest neighbour sites of the lattice. 

In this paper we examine the model in the range of 
parameters t,a « U . The value of a is usually smaller 
than that of t, however both these quantities are of 
the same order. Indeed, in systems described by the 
Hubbard- like models, it has been found that | a/t |« 0.3 
1^,^. In our studies we impose the following condition: 
—t<a<t (for a — 0, this model reduces to the ordinary 
FKM). 

General outline of the perturbation scheme. 

The perturbative scheme we use here can be applied 
to a general class of (lattice) Hamiltonians (defined on 
A C Z'') of the following form 



HA{t)=Ho^A + rVA, 



(4) 



where the unperturbed Hamiltonian i?o,A is a classical 
operator (i.e. it is diagonal in a basis being the prod- 
uct of bases on all lattice sites) with degenerate ground 
states. (Obviously the chFKM, introduced by the formu- 
las above, belong to this class). Our purpose is to exam- 
ine the effect of a quantum perturbation rVX, where r is 
a small parameter. In other words, we want to (block) 
diagonalize the Hamiltonian and find its ground states. 

To accomplish this task we are looking for an unitary 
transformation U{r), which (block) diagonalizes the full 
Hamiltonian. In most cases, finding out such a trans- 
formation exactly is a hopeless job. More constructive 
method is a perturbative treatment, consisted in "killing" 
the off-diagonal part of perturbation up to some finite 
power of the parameter r: 



H{r) = Ho+rV ^ H^"'^ (r) + r 



where Hq (r) is block-diagonal up to the order ri -|- 1 in 
r. Following this way one can determine an explicit for- 
mula for a diagonal part of the perturbed Hamiltonian 

r(") 



Results based on various perturbative schemes have 
been previously obtained for the simplest versions (i.e. 
without correlated hopping term) of the FKM and Hub- 
bard models: The method we applied in our 
studies was developed in [p^-p^ for Hamiltonians of the 
form (^. Their components can have quite general form; 
it is sufficient that both of them are sums of finite-range 
operators (or even infinite-range but exponentially de- 
creasing with distance). Moreover, it is assumed that 
Hq is expressible by translationally invariant m-potcntial 

SB- 

The technique developed in jT^-|T^ has also two other 
important aspects. First, since we restrict ourselves 
to the low-temperature region of the phase diagram, 
we need to diagonalize only a low- energy part of the 
Hamiltonian, what considerably simplifies the calcula- 
tions. Second, special care is taken to the form of the 
transformed Hamiltonian: it is formulated as a sum of lo- 
cal operators. It is necessary to obtain uniform estimates 
(i.e. independent of volume) and, as a consequence, to 
establish the convergence of the whole procedure, and (in 
some cases) to examine orderings emerging in the system. 

All this formalism, its background, achievements and 
limitations can be found in |16-19|. It must be stressed 



that many analogous results have been also obtained by 
Kotecky and coworkers . 

Effective Hamiltonian. The Hilbert space of the 
whole system Ti.A is a tensor product: Ti.A = ^^Ti-x- 

kSA 

Every Tix is spanned by the states: \wx,nx)- There are 
four base vectors:|0, 0), |1,0), |0, 1) and |1,1). The cor- 
responding energies are: 0; —^f, — /ie; U ~ ~ /ig. 

Let us begin our analysis from the classical part of the 
Hamiltonian. It is identical to a classical part of the Hub- 
bard model and the FKM, and it is well known The 
phase diagram consist of four regions. In region /, de- 
fined by fii < 0, /ie < 0, all sites are empty. In two twin 
regions Hi, He given by conditions: 7/^: /i^ > 0, /i^ > /ie, 
fie < U (for He, one should interchange the subscripts 
i and e) all sites are in the |1,0) (corresp. |0,1)) state. 
In the region ///, given by: fii > U,fj.e > U, all sites 
are doubly occupied. The most interesting situation is in 
the neighbourhood of the fMi = /ie line between regions 
Hi and He, which corresponds to the half-filled band, 
where there is a macroscopic degeneracy. We will anal- 
yse mainly this region. 

After some relatively straightforward but lengthy cal- 
culations, performed partially with an aid of symbolic 
computation programs, we have obtained for d — 2 the 
following effective Hamiltonian up to the second order of 
the perturbation theory, i.e. up to terms proportional to 
(7-3: 



Hq (r), called the effective Hamiltonian. 
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where: - the classical one-half spin on the lattice 
site i; it is related to the variable Wi by the formula: 
Si = {wi — l)/2; tci = t + a; Bz,ijk ~ "bent" triples of 
spins i, J, k (i.e. the angle between bonds ij and jk is tt/ 
2); S^^ijk - "straight" triples; V^^iju is a 2 x 2 plaquette 
on the lattice; h = fii — fie- 

Remark 1. In a general case, where some configura- 
tions of the system are not necessarily half- filled, there 
are also present projections onto half- filled states. For 
simplicity, we omit this aspect here. An expression for 
the general case we plan to present in an extended version 
of this paper. 

Remark 2. For a = 0, we should obtain an effec- 
tive Hamiltonian for the ordinary FKM. Comparing (|^) 
for a = with analogous expression in [|l9|, Table 2, 
we observed full consistence with exception of the con- 
stant term (we claim that authors have omitted term 

'Tu^^txyzw} (^'^ their terminology)). But this term is im- 
portant only for absolute values of energy; it neither af- 
fect the differences of energies, nor the phase diagram. 

Remark 3. The effective Hamiltonian (||) written in 
the spin variables corresponds to the Ising-like model 
with (dominating) antiferromagnetic interactions. How- 
ever there is an important new aspect, as compare with 
the simplest FKM, that comes from the correlated hop- 
ping term (i.e. for a ^ 0). This is a presence of terms 
with odd numbers of the spin operators. As a result the 
symmetry h —h (present in FK and Hubbard models) 
does not longer hold. An additional term proportional to 
the sum of spin variables plays a role of a suplementary 
external field. The other odd terms can be regarded as 
generalized fields. 

A physical explanation of the reason why these new 
terms with odd numbers of spin operators emerge is 
straightforward. Let us first focus on the linear term. 
Since in the chFKM the hopping rate depends on sites 
occupations, it is energetically favorable if all occupa- 
tion numbers have one, out of two possible values: 1 if 
a > (as hopping amplitude between two occupied sites 
is equal to < -|- a) and if a < 0. (Note that for the sim- 
plest FKM a = so the linear term does not contribute 



to the Hamiltonian, as it is proportional to a). The odd 
terms of higher order enter in a more subtle way, as they 
involve three or more neighbouring sites, but the princi- 
pal rule stay the same: it favours an occupation number 
equal to 1 or 0, for a > or a < 0, respectively. 

The phase diagram. Since there is no general 
method of finding ground-states for classical Hamiltoni- 
ans (as far as we know), we looked for the ground states 
of the Hamiltonian (||) by minimizing energy in some set 
of "trial" configurations (the method of restricted phase 
diagrams, [^,^ ) . We took a set of all periodic configura- 
tions, having elementary cells up to 12 sites (there are 
2000 such nonequivalent configurations), however it ap- 
peared that all configurations that emerged in the phase 
diagram, have no more than 5 sites per elementary cell. 
This observation led us to the conjection, that within the 
assumed perturbation order our results are exact, i.e. we 
claim that other configurations are absent. We hope that 
in a future we will be able to find a rigorous proof that 
our configurations are true minimizers. 

The phase diagram in variables (a/t, h) is displayed in 
Fig. 1. It can be noticed that for a predominant set 
of model parameters the sequence of phases agrees with 
that one found for the ordinary FKM (although the val- 
ues of h separating subsequent phases depend strongly 
on a). However, it is remarkable that for h < and 
a e]a_,a+[, where a+ = (-4- V2)/7 « -0.773459, 
a_ = (-4 + V^)/7^ -0.369398, a new type of order- 
ing (labelled by (4)) appears, instead of the three phases 
(3,5,6). Consequently, the diagram is clearly asymmetric 
with respect to the horizontal axis h — 0. 
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FIG. 1. Schematic phase diagram of the effective Hamilto- 
man ^ of the correlated hopping FKM. Phases represented 
by various arrangements of the ions are depicted in Fig. 2. 



Another interesting feature of the diagram is presence 
of the junction point for h — 0, a — —t, where all lines 
separating various phases for /i > 0, as well as most of 
the lines for h < 0, join together. It worthwhile to no- 
tice that the exact solution has been given just for this 
characteristic symmetry point [n5[. 
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FIG. 2. Configurations of the ions (marked by heavy dots 
•) corresponding to phases displayed in Fig. 1. Phases la- 
belled by numbers with prime (e.g. 1', 3', 4' etc.) have mirror 
configurations with respect to those without prim, i.e. lattice 
sites occupied by the ions are then interchanged with those of 
unoccupied by the ions. 

Summary. The results obtained in this paper one can 
summarize as follows. First, the effective Hamiltonian 
of the chFKM in the second order perturbation theory 
(i.e. up to terms proportional to U~^) has been found, 
and then its ground state phase diagram has been con- 
structed. It has become evident from our studies that the 
correlated hopping term modifies substantionally the ef- 
fective Hamiltonian and consequently the phase diagram 
of the simplest FKM. The main new feature of the ef- 
fective Hamiltonian is presence of odd parity terms. As 
the result the phase diagram becomes asymmetric with 
respect to change of a sign oi h. In particular, a new 
ordered phase, that does not exist for the simplest FKM, 
has been found for a certain interval of the parameter a. 

The possible directions for further studies that emerge 
from our results are: taking into account subsequent 
terms of perturbation theory, investigation the system at 
low, but nonzero temperatures and inclusion additional 
small terms to quantum part of the Hamiltonian (for in- 
stance, we will allow hopping of the ions with a small 
amplitude ti << t^, thus obtaining strongly asymmetric 
Hubbard model with correlated hopping). 
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